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Non-Gaussianity in the Cosmic Microwave Background Temperature Fluctuations 

from Cosmic (Super-) Strings 
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We compute analytically the small-scale temperature fluctuations of the cosmic microwave back- 
ground from cosmic (super-) strings and study the dependence on the string intercommuting prob- 
ability P. We develop an analytical model which describes the evolution of a string network and 
calculate the numbers of string segments and kinks in a horizon volume. Then we derive the prob- 
ability distribution function (pdf) which takes account of finite angular resolution of observation. 
The resultant pdf consists of a Gaussian part due to frequent scatterings by long string segments 
and a non-Gaussian tail due to close encounters with kinks. The dispersion of the Gaussian part is 
reasonably consistent with that obtained by numerical simulations by Fraisse et al.. On the other 
hand, the non-Gaussian tail contains two phenomenological parameters which are determined by 
comparison with the numerical results for P = 1. Extrapolating the pdf to the cases with P < 1, 
we predict that the non-Gaussian feature is suppressed for small P. 



I. INTRODUCTION 

The imprint of cosmic strings on the cosmic microwave 
background (CMB) has been widely studied. Although 
cosmic strings are excluded as a dominant source of the 
observed large-angular- scale anisotropics [![, they could 
still be observable at small scales [U, S HI with new ar- 
cminutc CMB experiments, such as the South Pole Tele- 
scope Q or the Atacama Cosmology Telescope 0. Be- 
cause the structure of a string network is highly non- 
linear, it would naturally induce a non-Gaussian fea- 
ture in the CMB fluctuations. In particular, a moving 
straight string produces discontinuities in CMB temper- 
ature, called the Kaiser-Stebbins effect pj, and a temper- 
ature gradient map has been suggested as a means for de- 
tecting such an effect § (see also ||). The non-Gaussian 
feature would also appear in the bispectrum [l(J, which 
has attracted much attention in the CMB community 
[Til ]. Fraisse et al. [|[ found that the probability distri- 
bution function (pdf) of the temperature fluctuations has 
a non-Gaussian tail and negative skewness. These non- 
Gaussian features may help us distinguish cosmic string 
signals from other secondary effects and hence enhance 
their observability 

Recently, cosmic superstrings have attracted much at- 
tention in the context of inflation in string theory fl2l.[l3|. 
Cosmic superstrings have properties different from con- 
ventional field-theoretic cosmic strings. One of the obser- 
vationally interesting differences is concerning the inter- 
commuting probability P. It can be significantly smaller 
than unity for superstrings while P = 1 is normally 
assumed for field-theoretic strings (but see [HI). This 
difference may be used to distinguish superstrings from 
field-theoretic strings observationally. 

In this paper, we compute analytically the pdf of the 
small-scale CMB temperature fluctuations and study its 
dependence on P. At small scales where the primary 
fluctuations are damped, only the integrated Sachs- Wolfe 



(ISW) effect is relevant. Because the contribution from 
loops was shown to be insignificant Q, we focus on the 
ISW effect of long string segments and kinks. We first 
present the basic formulae for the temperature fluctu- 
ations induced by long string segments and kinks (l5j 
(section HI)) , and follow the evolution of the number den- 
sities of segments and kinks by combi ning and extending 
a velocity-dependent one-scale model [l6l. [T7j and a kink 
model [18[ (section |III|) . Then, in section ITVl we derive 
the pdf showing that the results in Q can be interpreted 
with our simple model. Also the P dependence of the 
pdf is presented and the non-Gaussianity is predicted to 
be suppressed for small P. Finally we summarize our 
results in section IVl 



II. TEMPERATURE FLUCTUATIONS DUE TO 
COSMIC STRINGS 

First we summarize the basic formulae for the CMB 
temperature fluctuations due to cosmic strings, following 
We denote the position of a cosmic string by r(t, a) 
where t and a are the time and position on the string 
worldshect, respectively The equations of motion and 
constraints in a flat spacctime are given by, 
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where the dot and prime denote the derivatives with re- 
spect to t and a, respectively. Photons obtain or lose 
their energies due to the gravitational field of cosmic 
strings. The temperature fluctuation in the direction n, 
A = AT/T, due to a straight segment with the length £ 
is written as 



A seg (n) = 8 



vT 



za scg Gfi arctan -, (4) 
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where G is the Newton constant, fi is the string tension, 
6 is the impact parameter of a photon ray, v = \ f\ is the 
velocity of the segment and 
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is a factor which represents the configuration of the seg- 
ment and the direction of the line of sight and can be 
positive or negative. In the limit that the impact param- 
eter is much smaller than the segment length, 6 -C £, Eq. 
(U]) is reduced to the well known formula, 



A seg (n) = 4tt 



-a ses Gfj,. 



(6) 



In fact, Eq. (J3J can be well approximated by Eq. (O for 
5 < £ while it approaches zero for S > £. Therefore a 
segment with the length £ has an effective cross section 

~e. 

On the other hand, a kink can be modeled as a non- 
smooth junction of two straight segments with different 
directions, r' [15j . Then the temperature fluctuation 
with the impact parameter 5 is 

s 

Akink(") = -4G^a k mk log 6(X k mk - 5), (7) 

-^kink 

where Lkink is a distance between kinks. The step func- 
tion 9(Lkmk — S) represents the effect that the fluctuation 
becomes negligible far from the kink, and akink represents 
the kink configuration, 



CUkink = n- p, p 



Ckink + O 



ffkink-0 



(8) 



where Ukink is the position of the kink and p characterizes 
the change of the direction of the string at the kink and 
represents the kink amplitude. 



III. ANALYTIC MODEL OF COSMIC STRING 
NETWORK 

In this section, we develop an analytic model which 
describes the behavior of a cosmic string network. First, 
the average string length £ and the rms velocity u rms 
are calculated using a velocity-dependent one-scale model 
pH , [l7T | . Then, the number of kinks in a horizon volume 
is calculated with an approach similar to [IH (see also 

Him, ni). 

In the velocity-dependent one-scale model, a string 
network is assumed to consist of straight string segments 
with the average length £ and the rms velocity u rms - This 
scale £ is also assumed to characterize the interstring dis- 
tance, that is, p scg = /i/£ 2 , where p scg is the energy den- 
sity of string segments, respectively. In terms of £, the 
number of segments in a horizon volume is expressed as 



1 



£3 H 3 



7 



(9) 



where H is the Hubble parameter and we defined 7 = 

i/m)- 

The evolution of the network of segments is deter- 
mined by the cosmic expansion and the energy loss 
due to loop formation. A loop formation can occur 
through the intercommutation of two segments or the 
sclf-intercommutation of a single segment. The charac- 
teristic timescale for loop formation is ~ ^/Pv Tms . For 
a universe with the scale factor a(t) oc t@, the evolution 
equations for 7 and u rms are given by [HI 



7 at 2 



s 7 - Pv. 



2 

rms 1 



dt 



= (1 - v^ ms )H [fc(w rms )7 ~ 2v rri 



(10) 
(11) 



where c is a constant which represents the efficiency of the 
loop formation and k(v lms ) « (2v / 2/7r)(l — 8w 1 6 ms )/(l + 
8i>j ms ) is the momentum parameter |lq |. Hereafter we 
assume a matter-dominated universe and set j3 = 2/3. 

It is known that a string network approaches a "scal- 
ing" regime where the characteristic scale grows with 
the horizon size [24| . This means that 7 and v rms are 
asymptotically constant in time. Here we assume that 
the scaling behavior is already realized by the recombi- 
nation time. From (p~0|) and (fTTj) . we obtain the scaling 
values of 7 and u rms neglecting their time derivatives. For 
small cP they can be approximately given as, 
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We see that small P, which means the inefficient loop 
formation, leads to large 7 and hence large A^ sog . From 
Eq. (p~2|) . we have the dependence p scg oc P _1 , which 
is consistent with the result of numerical simulations in 
[23l | while [13] obtained a relatively weaker dependence 
on P. Actually there is no consensus on the dependence 
on P and we argue the effects of this ambiguity on the 
pdf later. 

Next, we consider the kink number evolution. Our 
approach is based on the idea of [H[ although the for- 
mulation is somewhat different. Small scale structure on 
strings including kinks is also considered in in a dif- 
ferent approach. Kinks are formed on string segments 
when they intercommute and, simultaneously, some of 
the existing kinks are removed through loop formation. 
Furthermore, kinks decay due to stretching by the cos- 
mic expansion and the emission of gravitational waves. 
Here we neglect the decay due to the gravitational wave 
emission and focus on the decay due to cosmic expan- 
sion since it is the most efficient decay process at the 
matter-dominated stage 

According to [25|, the kink amplitude, p = \p\ (see Eq. 
©), decays with cosmic expansion as p(t) = pf(£/tf)~ c , 
where if and pt arc the formation time and the amplitude 
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at the formation, respectively, and 
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We count the number of kinks with amplitude p m \ n < P < 
Pmax where p m in and p max are free parameters. Later we 
show that we need only the ratio Pmax/Pmin for our cal- 
culation and it will be determined by comparing our pdf 
for P = 1 with that obtained by the numerical simula- 
tions ||. Even if a kink is formed with pi > p m i n , the 
cosmic expansion reduces the amplitude gradually and 
eventually it is no longer counted as a kink after a time 
determined by p(t) — p m i n . Therefore, the kink number 
in a comoving volume V(t) oc a 3 (t) is given by the fol- 
lowing integral of the formation rate dN{ orm (t, p) / dtdp, 



Mcink = 



dp 



dt 



*o(p) 



dN lorm (t,p) 
dtdp 



(14) 



where to(p) = ^(p/Pmax) 1 ^, and a barred quantity is a 
value in the comoving volume V(t). 

The formation rate of kinks, which is assumed here to 
be independent of p, is proportional to the loop forma- 
tion rate, dN\ oov / dt. Because the loop formation rate 
determines the rate of the loss of the energy of string 
segments, wc have [1^ 



dNf mm (t) _ q dNi oop (t) _ qcPv rms V 



dtdp 



Pn 



dt 



(15) 



where q is a constant which represents the efficiency of 
the kink formation and a is the average loop length in 
units of £. Performing the integrations in (fT4|) . the kink 
number in a horizon volume iV k i nk = -^kink/(V"i? 3 ) is 



-^kink 



2qcPv rms j 4 e 



3a 



Pn 



Pn 



(16) 



where we have assumed eC 1. Because -/V k i nk is indepen- 
dent of time, the kink number is also scaling. This means 
that the average distance between kinks, L k j nk , evolves 
in proportion to the horizon scale. In fact L k ; n k is given 
by 



L 



kink 



Nscgt 



1 



A" 



iV ki „k KH 
iVkinkT Mcink 



N s , 



T 



(17) 
(18) 



where the normalized linear kink density, K, is constant 
in time. Thus we have expressed the numbers of string 
segments and kinks in a horizon volume as functions of 
P. 



IV. PDF OF CMB FLUCTUATIONS 

Based on the elementary processes presented in section 
im and the network evolution model in section IIII1 we 



calculate the pdf of the CMB temperature fluctuations 
due to string segments and kinks. 

A photon ray is scattered by segments many times 
through its way from the last scattering surface to an ob- 
server. Hence the temperature fluctuation would behave 
like a random walk and the pdf from segments would be 
approximated by the Gaussian distribution. If we treat 
a segment as a particle with the cross section £ 2 as we 
discussed below Eq. ([6]), the optical depth is 



3 e 2 



N seg H*£ 



dz 



H(l + z) 



7log(l + z rec ), (19) 



where z rec « 1100 is the redshift at the recombination. 
This is estimated as 7log(l + z rcc ) « 16 for P = 1 and 
larger for smaller P. Although the temperature change 
at each scattering is different depending on the factors 
a seg and v, it would be a good approximation to esti- 
mate the dispersion of the pdf using their statistical av- 
erages. Therefore, remembering Eq. ([S]), the dispersion 
is evaluated as, 



2tt- 



2 

e 



27ra seg \/l0g (1 + Zrec) 



7Tv/2 

3c\P 



1/4 



G», (20) 



where we have set v = \ 
in the third equality, and A sog and a seg should be un 



and substituted Eq. |T 



derstood as their statistical averages, 



< A 2 0g > and 

< a 2 eg >, respectively. Here it should be noted that 

the PDF from string segments should be, strictly speak- 
ing, the binomial distribution. However, for the number 
of trials evaluated above (~ 16), the deviation of the 
binomial distribution from the Gaussian distribution is 
negligibly small. 

Next, let us consider the contribution from kinks. The 
temperature fluctuation depends on the impact parame- 
ter as given by ([7]). Solving for S as a function of A, we 
have, 



5(A) = L kink e-I A I/ 2A <\ 
A = 2a kink G^. 



(21) 
(22) 



Therefore the differential cross section with the temper- 
ature fluctuation A can be written as 



dcr k i, 



dA 



dA 



5 2 (A) 



T 2 

_ ^kink -IAI/Aq 

" A 



(23) 



where a k i nk should again be understood as its statistical 
average of the kink configuration and cr kink should not 
be confused with the coordinate on a string in section HT1 
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dPki 



kink 



dA 



.do- 



kink 



dz 



dA H(l + z) 



(24) 



Then the pdf of temperature fluctuations due to kinks is Here we just treat them as phenomcnological parameters 

and put akink = 4.5 and p max /p m i n = 2.3 to make (p2"9"]) 
and (|30[) consistent. 

Before we discuss the P dependence of the pdf, let us 
consider the effect of finite angular resolution of CMB 
observation. As we saw above, a ray has to pass nearby 
a kink to have a large temperature fluctuation but it 
may not be resolved if the angular resolution is finite. 
This effect can be taken into account by assuming that 
the impact parameter 5 cannot be smaller than a certain 
value S m i n (z, 8) determined by the redshift of a kink and 
(3/2) v / 3v / 2/7rcP the angular resolution 8, 

(25) , , 

<WM) = 8d A (z) = 28H- 1 (z)(VT+^-l), (31) 



= ^ e - |A|/Ao log(l + z roc ) 

The normalization factor can be evaluated as 
7 2 log (1 + z rec ) 



.4 



KA 

37ralog (1 + z r 
8qa kink Gn 



=) ( 3V2 

1 TTcP 



3/2 



Pn 



Pu 



where we have used Eqs. (|12|) . ([16]) and (fl8[> in the second 
equality Thus we have a pdf of the form, 



dPtot 


dP G dP NG 


(26) 


dA 


dA dA ' 


dP G 


1 e -A 2 /2 CT 2 


(27) 


dA 


v/2^cr 


dP NG 
dA 


Af HA|/A 0j 


(28) 



where a, Ao and A are given by Eqs. PO]) , (|2"2"]) and 
([25]) . respectively. dP G /dA is the Gaussian part due to 
frequent scatterings by string segments, and dP^o/dA is 
the non-Gaussian tail due to rare scatterings by kinks. 
Here, because dP^c/dA < 1 as we see just below, we 
have normalized dPo/dA as dA dPt^/dA = 1. 
In the limit P — > 1, we have 



A = 2a kink G^, 



Pr 



Pn 



-5.1 



(Gm)" 1 , 

(29) 

0.23 and a = 0.1 as their 



Here we used aw = 1/ v2 for the 



where we have set 

standard values [2^. ucn wo uo™ ct SO g 
statistical average assuming a random distribution of r ' 
and r. Contrastingly, the statistical average of c^ink can- 
not easily be obtained because kinks evolve in time and 
their distribution is nontrivial. This problem is closely 
related to the number count of kinks discussed in section 
Mil and we will postpone it to the future work. 

On the other hand, the pdf from numerical simulations 
can be also described as Eqs. ((26]), ([27]) and J2SJ with 



12Gn, A sim « 0.03(G A1 )~ 1 , A , ; 



9G[i. (30) 



First, we note that the dispersion of the Gaussian part 
is well reproduced without any adjustable parameters. 
This would imply that our interpretation of the Gaus- 
sian part as frequent scatterings by segments is reason- 
able. Next, to compare the non-Gaussian part of Eqs. 
([29]) with ([30]) . we must specify the values of akink and 
Pmax/Pmin- To estimate these parameters from the first 
principles is, however, beyond the scope of the present pa- 
per. This will be discussed in a forthcoming paper [271 ] . 



where d^(z) is the angular diameter distance. The 
largest fluctuation which can be generated by a kink at 
z is, then, 



2A log 



^kink(z) 
S m in(z,8 



= -2A log[26K(y/T+z-l)]. (32) 

Solving this in terms of z, we obtain the largest redshift 
of kinks which contributes to a specific value of A, 



c(A, 



3 -|A|/2A 

28K 



1, Zr, 



(33) 



where this is bounded by z roc because we are considering 
the ISW effect. Then the pdf modified by a finite angular 
resolution is obtained by changing the integration range 
of Eq. JH from [0,^ rcc ] to [O,w(A,0)], 



jprcs 
"-'kink 

dA 



k(A,S) 



-/Vkink-ff 



, da 



kink 



dz 



dA H(l + z) 



?T_ e -|A|/Ao 



KA 
log/l 



,-|A|/2A 

28K 



1 



,-|A|/2A 



2e -|A|/2A +2K8 



(34) 

where we have assumed z max (A, 8) < z roc . In the limit 
of an infinite resolution, this reduces to Eq. (|24| . and in 
the opposite limit, we have 



dPki: 



37 2 



-3|A|/2A 



dA 28K 2 A 



(35) 



The effect of a finite resolution is important for |A| > 
Ai where Ai is defined by e~ Al/2A ° /(28K) = 1, that 
is, Ai = -2A log(2^). Setting 8 = 0.42' which was 
adopted in Ai is estimated as w 68G/i for P = 1 
and larger for smaller P. Thus the non-Gaussian tail 
steepens slightly for large |A| but the pdf of 0] is still 
well reproduced. 

In Fig. [TJ the dependence of the pdf on the intercom- 
muting probability P are shown with angular resolution 
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FIG. 1: Dependence of the pdf with angular resolution 9 = 
0.42' on the intercommuting probability P (thick lines). The 
respective Gaussian parts are plotted with thin lines for com- 
parison. For P = 1 and 10 -0 ' 25 , the pdfs deviate significantly 
from the Gaussian distribution. For P < 10~ 0,5 , pdfs are 
almost Gaussian. 




A/Gn 

FIG. 2: Dependence of the pdf with P — 1 on angular reso- 
lution. The Gaussian part is also plotted for comparison. 



6 = 0.42'. The respective Gaussian parts are also plot- 
ted. As we see, as P decreases, the Gaussian dispersion 
increases and the contribution of the non-Gaussian tail 
is suppressed. For the case of P = 1, the pdf is almost 
Gaussian for |A| < 50G/i while non-Gaussian tails can 
be seen for larger |A|. Contrastingly, the pdfs are al- 
most Gaussian for a wide range of |A| in the cases with 
P < 10~ " 5 . Thus the non-Gaussianity could be a probe 
of the cosmic string property, P. 

Fig. [5] shows the dependence on angular resolution 
with P = 1. We see that for a typical angular resolu- 
tion, 5', of future observations such as Planck [28], non- 
Gaussian feature is highly suppressed even for P = 1 
because kinks can not be resolved. Thus we would need 
observations with an arcminute resolution. 

Note that our pdf is symmetric for positive and nega- 



tive A and cannot reproduce the non-zero skewness re- 
ported in This is because we have assumed the long 
segments to be straight. In fact, for a straight segment, 
temperature fluctuation is symmetric between positive 
and negetive values as is seen from Eqs. ([4]) and ([5]). In 
this case, even if the number of scatterings is relatively 
small and the Gaussian approximation is not valid, skew- 
ness does not appear. 

However, skewness would appear if we take the cur- 
vature of se gme nts and its correlation with velocity into 
account [271 129||. The reason is as follows. 

First, if a segment has a curvature, the symmetry of 
temperature fluctuation between positive and negative 
values is broken. The asymmetry depends on the angle 
between the velocity and curvature vectors. Therefore 
the curvature cannot induce a skewness itself because the 
asymmetry would be canceled out by multiple scatterings 
of strings with various configuration of velocity and cur- 
vature vectors. However, if there is a correlation between 
velocity and curvature, there would be a nonzero expec- 
tation value for the asymmetry. This is the mechanism 
we believe the skewness is induced from. 

To estimate the skewness, we need to extend the fluc- 
tuation formula Eq. (g]) taking the curvature of a seg- 
ment into account. With the nonzero expectation value 
of asymmetry, the deviation from the Gaussian distribu- 
tion would lead to a nonzero skewness. The extcntion 
of Eq. ((3]) and the evaluation of the deviation from the 
Gaussian distribution due to the finite number of scat- 
terings will be discussed in a separate article (2?| ■ Never- 
theless, it would be surprising that most of the features 
of the pdf obtained by the numerical simulations Q can 
be interpreted by our simple model with just straight 
segments and kinks. 



V. DISCUSSION AND SUMMARY 

In this paper, we have computed analytically the pdf 
of small-scale CMB temperature fluctuations due to cos- 
mic (super-)strings with a simple model with straight 
segments and kinks. Our purposes were to interpret the 
results of numerical simulations in Q and study the ef- 
fect of the string intercommuting probability P. We have 
combined and extended a velocity-dependent one-scale 
model and a kink model to calculate the numbers of 
string segments and kinks in a horizon volume consis- 
tently Thus obtained pdf consists of a Gaussian com- 
ponent due to frequent scatterings by string segments 
and a non-Gaussian tail due to close encounters with 
kinks. The dispersion of the Gaussian part obtained 
by numerical simulations [3[ is well reproduced without 
any adjustable parameters. On the other hand, the non- 
Gaussian tail contains two phenomenological parameters 
and we determined them by comparing it with that of 
the numerical result for P = 1 by Fraisse et al. Q. Then 
we clarified the P dependence of the pdf and found that 
the non-Gaussian tail diminishes as P decreases. 
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Let us argue the ambiguity in our string network model 
in section IIII1 The evolution of large- and small-scale 
structure has not been well understood either analyti- 
cally or numerically. In particular, the dependence of 7, 
i> rms and K on P is quite important to derive the P de- 
pendence of the pdf by our formalism. As we pointed 
out below Eq. (|12p , the dependence of 7 in our network 
model is consistent with that of [23| while [l7| claims a 
relatively weaker dependence. However, they are consis- 
tent in that a small P results in large 7, N seg and p se g- 
Then, from the second equation of Eq. (f2Tj)) . it would be 
robust that the dispersion of the Gaussian part increases 
as P decreases suppressing the non-Gaussian tail, assum- 
ing that v lms would not differ significantly from 1 / \[2 (sec 
Eq. (O). Anyway, the parameters of the pdf can be cal- 
culated by our formalism once 7, i> rms and K are given as 
functions of P. Thus it would be interesting to calculate 
the pdf using those functions obtained by other models 
and numerical simulations. 

Our pdf is contributed only from the ISW effect of cos- 
mic strings. Although the primary temperature fluctu- 
ations are substantially damped at small scales we con- 
sider here (~ 0(1) min), other secondary fluctuations 
such the as Sunyaev-Zel'dovich effect would become im- 
portant depending on the value of G[i. To discuss fur- 
ther observational prospects of cosmic (super-)strings, it 
would be necessary to compare contributions from vari- 



ous secondary fluctuations. 

Also it is important to compute other observational 
quantities with our simple formalism. In particular, the 
power spectrum contributed from cosmic strings has been 
calculated by many authors @, H, [1] and the comparison 
with them would further allow us to check the applica- 
bility of our formalism. As to non-Gaussianity, observa- 
tionally more interesting quantities than the pdf would 
be higher-order correlation functions. The work along 
this direction is in progress [27| . 
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